How long can a quantum memory withstand depolarizing noise? 
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We investigate the possibilities and limitations of passive Hamiltonian protection of a quantum 
memory against depolarizing noise. Without protection, the lifetime of a single qubit is independent 
of N, the number of qubits composing the memory. In the presence of a protecting Hamiltonian, 
the lifetime increases at most logarithmically with N. We construct an explicit time-independent 
Hamiltonian which saturates this bound, exploiting the noise itself to achieve the protection. 
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A cornerstone for the majority of applications in quan- 
tum information processing is the ability to reliably store 
quantum information, protecting it from the adversarial 
effects of the environment. Quantum Error Correcting 
Codes (QECC) achieve this task by using a redundant 
encoding and regular measurements which allow for the 
detection, and subsequent correction, of errors 0, 0, 0, H| • 
An alternative approach uses so-called protecting Hamil- 
tonians [H, 0], which permanently act on the quantum 
memory and immunize it against small perturbations. Its 
most attractive feature is that, in contrast to QECC, it 
does not require any active action on the quantum mem- 
ory, just encoding and decoding operations at the time of 
storing and retrieving the information. Whereas this ap- 
proach may tolerate certain types of perturbation 0, 0] , 
it is not clear if it is suitable in the presence of depolar- 
izing noise, something which QECC can deal with. 



In this Letter, we give a complete answer to this ques- 
tion. More specifically, we consider the situation where 
a logical qubit is encoded in a set of N qubits and al- 
lowed to evolve in the presence of depolarizing noise and 
a protecting Hamiltonian. The goal is to find the strategy 
delivering the longest lifetime, r, after which we can ap- 
ply a decoding operation and reliably retrieve the original 
state of the qubit. By adapting ideas taken from [9(, it 
is established that the lifetime cannot exceed log TV. An 
analysis of the case in which no protecting Hamiltonian is 
used presents markedly different behavior depending on 
whether we intend to store classical or quantum informa- 
tion. Finally, we construct a static protecting Hamilto- 
nian that saturates the upper bound t ~ O(logiV). To 
this end, we first show how to achieve this bound using a 
time-dependent Hamiltonian protection which emulates 
QECC. We then introduce a clock gadget which exploits 
the noise to measure time-similar to radiocarbon dating- 
thus allowing us to simulate the previous time dependent 
protection without explicit reference to time. 



We consider a system of N qubits, each of which is 
independently subject to depolarizing noise at a rate r. 



The total state evolves as 



p(t) = -i[H(t),p(t)]-r 



n=l 



(1) 

Furthermore, we shall allow for an arbitrary encoding of 
the initial state as well as a final decoding procedure to 
recover the information. 

Protection limitations. — Using purely Hamiltonian 
protection, a survival time of r ~ 0(log N) is the maxi- 
mum achievable. Intuitively, this is due to the fact that 
the depolarizing noise adds entropy to the system at a 
constant rate, while any reversible operation (i.e., Hamil- 
tonian/unitary evolution) will never be able to remove 
this entropy from the system. Rather, in the best case, it 
can concentrate all the entropy in a subsystem, keeping 
the remaining part as clean as possible. This entropic 
argument was first presented in where the authors 
investigated the power of reversible computation (both 
classical and quantum) subject to noise in the absence of 
fresh ancillas. To this end, they considered the informa- 
tion content I(p) = N — S(p) of the system, with N the 
number of qubits and S(p) = — tr(plog 2 p). The infor- 
mation content upper bounds the number of classical bits 
extractable from p, and thus ultimately also the number 
of qubits stored in p. While the original statement about 
the decrease of I(p) is for discrete-time evolution, it can 
be straightforwardly generalized to the continuous time 
setting of Eq. ^ , where it states that 



dt 



< -rl(p) , 



which yields that the information of the system is smaller 
than i after a time . 

2 r 

Having established an upper bound for the scaling of r 
with N, let us analyze whether this bound can be reached 
under different circumstances. We start out with the 
simplest case where we use no Hamiltonian protection 
(i.e., H = 0) and show that r is independent of N; that 
is, no quantum memory effect can be achieved. For that, 
we note that the effect of Eq. |lj on each qubit may be 
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expressed in terms of a depolarizing channel 



£ t ( P ) = X(t)p- 



1 



A(t))| 



the 



1 

This re 



where A(i) = e _rt . For i > i c i, where 
resulting channel is entanglement breaking [lC 
mains true if one incorporates encoding and decoding 
steps. It is simple to prove that for entanglement break- 
ing qubit channels, the average fidelity [ll| is upper- 
bounded by 2/3. Thus, the lifetime r is smaller than 
t c \ = ln3/r, which is independent of N. 

The previous argument docs not apply to classical in- 
formation, for which an optimal storage time that is log- 
arithmic in N may be achieved. The classical version 
of Eq. [TJ taking H{t) = 0, is a system of N classical 
bits subject to bit flipping noise (each bit is flipped at 
a rate r/2). In this case, encoding in a repetition code, 
and decoding via majority voting, yields an asymptoti- 
cally optimal information survival time O(logiV). Using 
optimal estimation [12| and this classical protocol in the 
encoding phase, the bound 2/3 may be asymptotically 
reached for the quantum case. 

Time dependent protection. — We will now use the 
ideas of QECC to build a simple circuit based model 
that reaches the upper bound on the protection time. 
This model assumes that unitary operations can be per- 
formed instantaneously, which is equivalent to having a 
time-dependent protecting Hamiltonian with unbounded 
strength; we will show how to remove both requirements 
later on. Instead of using a repetition code, we encode the 
qubit to be protected in an I level concatenated QECC 
(he., I levels of the QECC nested into each other), 
which requires N = d l qubits, where d is the number of 
qubits used by the code. Each level of the QECC can pro- 
vide protection for a constant time t plo t < t c \, and thus, 
after t plot one layer of decoding needs to be executed. 
Each decoding consists of a unitary t/dec on each <i-tuple 
of qubits in the current encoding level; after the decoding, 
only one of each of the d qubits is used further (Fig. [1]). 
The total time that such a concatenated QECC can pro- 
tect a qubit is given by t pmt l = t prot \og d N ~ O(logiV), 
as in the classical case. 

Time-independent protection. — In the following, we 
show that the same logiV protection time which we 
can achieve using a time-dependent protection circuit 
can also be obtained from a time-independent protect- 
ing Hamiltonian. The basic idea of our construction is 
to implement the time-dependent Hamiltonian presented 
before in a time independent way. To this end, a clock 
is built which serves as control. The time-independent 
version performs the decoding gates conditioned on the 
time estimate provided by the clock. In order to obtain a 
clock from ([I]) with a time-independent H, we will make 
use of the noise acting on the system: we add a number, 
K, of "clock qubits" which we initialize to and let 

the depolarizing noise act on them. The behavior of the 




FIG. 1: Decoding a nested QECC. The "discarded" qubits 
carry most of the entropy and are not used further. 

clock qubits is thus purely classical; they act as K classi- 
cal bits initialized to 1 which are being flipped at a rate 
r/2. Thus, the polarization k, defined by the number of 
"1" bits minus the number of "0" bits has an average 
expected value of k(t) = Ke~ rt at time t. Conversely, 
this provides the time estimate 



t(k) = min 



f]n{K/k) 



(2) 



Here, i max is the maximum time for which we expect the 
estimate to be reliable which depends on K and the pre- 
cision of the estimate, cf. §5§ below. Particular realiza- 
tions of this random process of bit flips can be described 
by a polarization trajectory k(t). Good trajectories are 
defined to be those such that 



|jfe(f) - k(t)\ < K^ 2+€ 



(3) 



for all < t < i max . For appropriate parameters £ max 
and < e < i, the following theorem states that almost 
all trajectories are good and can provide accurate time 
estimates. 

Theorem (Depolarizing clock): For K > 16, good 
trajectories have a probability 



P [k(t) good traj.] > 1 - K 



ft,, 



exp[-3if 2e /8] 



cxp[K 2e 



(4) 



Furthermore, for any good trajectory k{t), the time esti- 
mate t returned by the clock will differ from the real time 
t by at most 



1 



r Kl/2-e 



> \t(Ht)) - t\ 



(5) 



(For fixed 5, this implies that t m ax will scale logarithmi- 
cally with the number of qubits.) 

Note that the theorem does not simply state that any 
time evolution will be outside ([3]) for an exponentially 
small amount of time (which is easier to prove) , but that 
there is only an exponentially small number of cases in 
which © is violated at all. Although the former state- 
ment would in principle suffice to use the clock in our 
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FIG. 2: (color online) A step-like trajectory in green illus- 
trates the two ways of leaving region ((3} of good trajectories 
(dashed lines): either a spin flip can take the polarization out 
of the marked region (thick blue), or polarization may leave 
region ([3| as time passes without a spin flip (red dots). 



construction, the stronger version of the theorem makes 
the application of the clock, and in particular the error 
analysis, more transparent and will hopefully lead to fur- 
ther applications of the clock gadget. 

Proof. To prove the theorem, note that each of the bits 
undergoes an independent exponential decay, so that the 
total polarization is the sum of K identical independent 
random variables. We can thus use Hocffding's inequality 
[3] to bound the probability of finding a polarization far 
from the expected average value k(t), 



Pr 



\k(t)-k(t)\ > K 1 ^ 



< 2e~ 



(6) 



This already implies that most of the trajectories vio- 
late ([3]) for no more than an exponentially small amount 
of time. To see why ((6j) implies that most trajectories 
are good trajectories, we bound the average number of 
times a trajectory leaves the region (j3]) of good trajec- 
tories. Since a non-good trajectory must leave ([3|) at 
least once, it is also an upper bound on the probability 
of non-good trajectories. Hence, it suffices to consider 
the average rate R(t) at which processes leave (|3|), and 
integrate over t to obtain a bound on the probability of 
trajectories which are not good. 

The rate at which a process leaves the set of good tra- 
jectories has two sources, as illustrated in Fig. [5] First, 
the system can undergo a spin flip, thus leaving the re- 
gion defined by J3]) vertically (rate R v ), and second, it 
can leave it horizontally if the time t passes the maxi- 
mum time allowed by (J3j) for the current value k(t) of 
the polarization (rate Rh)- A vertical leave can occur 
only if \k(t) - k(t)\ > K l ' 2+e - 2 > K 1 / 2+e /2, provided 
K > 16 (a spin flip changes k(t) by ±2). Eqn. ([6|) thus 
gives an average bound 

R v (t) < Kre- K2e / 8 . 

A horizontal leave can only occur at discrete times ex- 



tremizing (|3|), 

t G T = {t : k(t) + K^ 2+e G N} , 

and the probability of a trajectory fulfilling k(t) = k(t) ■ 
K 1 / 2+e may again be bounded using ((6]), such that 



Rh(t) < 2e 



-K 2 '/2 



The inequality (j4]) follows immediately by integrating 
R h (t) + R v {t) from to t max . 

Assuming that k(t) corresponds to a good trajectory, 
the accuracy of the time estimate ((2]) may be bounded 
by applying the mean value theorem to k: 



!*(*(*)) - *l 



\k(i(k(t)))-k(t)\ 

\k (^interm)| 



< 



K' 



K 



□ 

Clock dependent Hamiltonian. — Let us now show how 
the decoding circuit can be implemented using the clock 
gadget. The circuit under consideration consists of the 
decoding unitaries (decoding the Ai'th encoded qubit 
in level I, acting on d qubits each); after a time interval 
(the time one level of the code can protect the qubit 



prot 



sufficiently well), we perform all unitaries U l ^ c at the 
current level I — note that they act on distinct qubits and 
thus commute. Each of these unitaries can be realized by 



T l,k 



for a time t = t^c 



applying a d-qubit Hamiltonian H' dcc 
Thus, we have to switch on all the H^ e for t E [U,ti 
tdec], where i; = lt pmt + {I - l)t dcc . 

In order to control the Hamiltonian from the noisy 
clock, we define clock times fc/ i0n = |_&(ti)J an d &/,off = 
\k(ti +idcc)l, and introduce a time-independent Hamil- 
tonian which turns on the decoding Hamiltonian for level 
I between k G [^ j0n , fcj.off], 



I 



Li 



H 



l,d L 



n. 



(7) 



The left part of the tensor product acts on the N code 
qubits, the right part on the K clock (qu)bits, and 

fez, off 

fe=fej, on W:c = (k+N)/2 

where x is an iV-bit string with Hamming weight w x . The 
initial state of the system is, as for the circuit construc- 
tion, the product of the encoded qubit in an Z-level con- 
catenated code and the maximally polarized state 
on the clock gadget. 

Error analysis. — We now perform the error analysis 
for the protecting Hamiltonian (JT]). In order to protect 
the quantum information, we will require that the er- 
ror probability per qubit in use is bounded by the same 
threshold p* after each decoding step is completed (i.e. 
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at t = ti + idee + f )■ We will restrict to the space of good 
trajectories, since we know from the clock theorem that 
this accounts for all but an exponentially small fraction, 
which can be incorporated into the final error probability. 

We will choose K large enough to ensure that the er- 
ror | > \t — t\ in the clock time satisfies 5 <C i P rot,idoc- 
In this way, we ensure that the decoding operations are 
performed in the right order [l6| and with sufficient preci- 
sion. We may thus account for the following error sources 
between ti + idec + 5/2 and ti+i + tdoc + 5/2: 

i) Inherited errors from the previous rounds which 
could not be corrected for. By assumption, these errors 
are bounded by pinhcr < P* ■ 

ii) Errors from the depolarizing noise during the free 
evolution of the system. The system is sure to evolve 
freely for a time i pro t — 5, i.e., the noise per qubit is 
bounded by p cvo \ < 1 - cxp[-r(t prot - 6)] < r(t prot - 5). 

iii) Errors during the decoding. These errors affect the 
decoded rather than the encoded system and stem from 
two sources: On the one hand, the time the Hamilto- 
nian is active has an uncertainty idoc ± 5, which gives 
an error in the implemented unitary of not more than 
exp[<5||-ff d ' ||] — 1. On the other hand, depolarizing noise 
can act during the decoding for at most a time tdec + 5. 
In the worst case, noise on any of the code qubits during 
decoding will destroy the decoded qubit, giving an error 
bound d(l — exp[— r(idcc + 5)]) < dr(td ec + 5). Thus, the 
error on the decoded qubit is 

Pdoc < exp[||F d ; c ||5] - 1 + dr(t dc c + 5) . 

Since the noise is Markovian (i.e. memoryless), the clock 
does not correlate its errors in time. In summary, the 
error after one round of decoding is at most B{p m \ lcr + 
Pcvoi)+Pdoc, which we require to be bounded by p* again. 
Here, B{p) is a property of the code, and returns the error 
probability of the decoded qubit, given a probability p of 
error on each of the original qubits; for example, for the 
5-qubit perfect QECC [3, B{p) < 10p 2 . 

We will now show that it is possible to fulfill the re- 
quired conditions by appropriately defining the control 
parameters. First, we choose p* < 1/40 to have the 



QECC [l4( work well below threshold. We may take 



iprot ~ and tdec := 43- • To minimize imprecision in 
the implemented unitaries, the decoding Hamiltonians 
are chosen of minimal possible strength, ||-ff dl l c || < 7—- 



Finally we take 5 



8tt 



. Inserting the proposed values 



in the derived bounds, it is straightforward to show that 

5(pi n her + Povol) + Pdcc < P* ■ 

The number of code qubits required is N := d l , with 
I '■= It — t +«d 1 ' ^ e rcc l u i re d logarithmic clock lifetime 
tmax = T and the precision 5 are obtained by taking e = 
1/6 and K := (^fj-) 3 , by virtue of Eq. © of the clock 
theorem. For any fixed r and p* , this allows a lifetime 
t - 0(log(N + K)). 

Conclusions. — In this paper, we have considered the 



ability of a Hamiltonian to protect quantum informa- 
tion from decoherence. While without a Hamiltonian, 
quantum information is destroyed in constant time, the 
presence of time-dependent control engenders protection 
for logarithmic time, which is optimal. As we have 
shown, the same level of protection can be attained with 
a time-independent Hamiltonian. The construction in- 
troduced a noise-driven clock which allows a time de- 
pendent Hamiltonian to be implemented without explicit 
reference to time. 

Since depolarizing noise is a limiting case of local noise 
models, it is expected that the time-independent Hamil- 
tonian developed here can be tuned to give the same 
degree of protection against weaker local noise models, 
although these models may admit superior strategies. 
For instance, noise of certain forms (such as dephas- 
ing) allows for storage of ancillas, potentially yielding 
a linear survival time by error correcting without decod- 
ing. In the case of amplitude damping noise, the noise 
itself distills ancillas so that the circuit can implement 
a full fault-tolerant scheme, which gives an exponential 
survival time, assuming that one can redesign the clock 
gadget to also benefit from these properties. 

Whether the same degree of protection can be ob- 
tained from a Hamiltonian which is local on a 2D or 3D 
lattice geometry remains an open question (l7j . How- 
ever, intuition suggests this might be impossible; the cru- 
cial point in reversibly protecting quantum information 
from depolarizing noise is to concentrate the entropy in 
one part of the system. Since the speed of information 
(and thus entrop y) t ransport is constant due to the Lieb- 
Robinson bound [la ] , the rate at which entropy can be re- 
moved from a given volume is proportional to its surface 
area, while the entropy increase goes as the volume. It 
thus seems impossible to remove the entropy sufficiently 
quickly, although this argument is not fully rigorous, and 
the question warrants further investigation. 
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